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Abstract. The long standing problem of a non-perturbative renormalization of a gauge field theoretical 
Hamiltonian is addressed and explicitly carried out within an (effective) light-cone Hamiltonian approach 
to QCD. The procedure is in line with the conventional ideas: The Hamiltonian is first regulated by suitable 
cut-off functions, and subsequently renormalized by suitable counter terms to make it cut-off independent. 
Emphasized is the considerable freedom in the cut-off function which eventually can modify the Coulomb 
potential of two charges at sufficiently small distances. The approach provides new physical insight into 
the nature of gauge theory and the potential energy of QCD and QED at short distance. 
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1 Introduction 

One of the most interesting subjects of particle and nu- 
clear physics is to understand the internal structure of 
hadrons. A powerful description of the hadronic structure 
is obtained through the light-cone wave functions. They 
encode all possible quark and gluon momenta, and helicity 
and flavor correlations in the hadron and are constructed 
from the QCD light-cone Hamiltonian pQ: 



H 



LC 



p+p- 



P 



where P = P° ± P z . The wave function tf'f, of a hadron 
b with mass satisfies the relation pQ: 



Ml 



The light-cone wave functions are expanded in terms of a 
complete basis of Fock states having increasing complex- 
ity. In this way the hadron presents itself as an ensemble of 
coherent states containing various numbers of quark and 
gluon quanta. For example, the positive pion has the Fock 
expansion: 

= $>|7r+)|n) = *lf /w ( Xi ,k ±i )\ud) 

n 

+ <2 / >i,kxi)M ? ) + ... , 

representing an expansion of the exact QCD eigenstate at 
scale A in terms of non- interacting quarks and gluons. The 
i n particles in a Fock state (n) have longitudinal light-cone 
momentum fractions Xi and relative transverse momenta 
k_u, with 



pO + pz 



i=l 



= 1 



The form of ^ n /jj(xi, kj_j) is invariant under longitudinal 
and transverse boosts; i.e., the light-cone wave functions 
expressed in the relative coordinates Xi and fcj_i are inde- 
pendent of the total momentum (P + , Pj_) of the hadron. 
The first term in the expansion is referred to as the va- 
lence Fock state, as it relates to the hadronic description 
in the constituent quark model. The higher terms are re- 
lated to the sea components of the hadronic structure. It 
has been shown that the remining terms of the light-cone 
wave function can be determined once the valence Fock 
state is known 2,3 . Explicit expressions are given in [3J- 
When starting in 1984 with Discretized Light-Cone 
Quantization (DLCQ) and with a revival of Dirac's 
Hamiltonian front form dynamics jS], all challenges of a 
gauge field Hamiltonian theory were essentially open ques- 
tions, particularly 

• the non-perturbative bound state problem, 

• the many-body aspects of gauge field theory, 

• regularization, 

• renormalization, 

• confinement, 

• chirality, 

• vacuum structure and condensates, 

just to name a few. The step from the gauge field QCD La- 
grangian down to a non-relativistic Schrodinger equation 
was completely mysterious. Now we know better £Q. 

We have understood, for example, that the chiral phase 
transition, in which the quarks are supposed to get their 
mass, is not the major challenge. The challenge is to un- 
derstand what happens after the phase transition, at zero 
temperature. 

The challenge is to understand the spectra of physical 
hadrons and to get the corresponding eigenfunctions, the 
light-cone wave functions. 
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We have learned how to partition the problem and how 
to shape our thinking in four major steps: 

£qcd 







M 2 \<P) 


(1) 














M 2 \*qq) 


(2) 




a- 






H e & \Vqq) 




E Wqq) 


(3) 



a- 

(ik + y(x) ) ^ (x) = E ^ h (4) 

The step from the Lagrangean £qcd to the light-cone 
Hamiltonian H^c had been given in pQ; all its matrix ele- 
ments are known and tabulated. Its reduction to an effec- 
tive light-cone Hamiltonian i?LC has been given in 3,6 
and will be sketched once again in Sec. 01 Some time will 
be spent in Sec. to discuss the single particle Hamilto- 
nian H e f{ in momentum and in configuration space. 

But first we must deal with the key issue of any gauge 
theory. We must say something about how to treat the 
logarithmic divergencies of the theory. We shall do it by 
conventional regularization and renormalization in Sees. El 
and[7| In the perturbative context of scattering theory the 
problem of regularization and renormalization was solved. 
In the non-perturbative context of a Hamiltonian it has 
not been solved. It is addressed in the next two sections. 

Preliminary and condensed versions of this were given 
in conference proceedings [7] . Several typos and misprints 
are removed, also including some mistakes. 

2 Regularization of the total Hamiltonian 

Canonical field theory with the conventional QCD La- 
grangian allows to derive the components of the total 
canonical four- momentum P+ Its front form version pQ 
rests on two assumptions, the light cone gauge A + = 
|S] and the suppression of all zero modes HE|. The front 
form vacuum is then trivial. 

I find it helpful to discuss the problem in terms of 
DLCQ PI3] . In the back of my mind I visualize an explicit 
finite dimensional matrix representation of the Light-Cone 
Hamiltonian as it occurs for finite harmonic resolution. 
Such one is schematically displayed in Fig. 2 of pQ. All of 
its matrix elements are finite for any finite x and k±. 

The problem arises for ever increasing harmonic reso- 
lution, on the way to the continuum limit: The numerical 
eigenvalues are numerically unstable and diverge logarith- 
mically jlOlllj . in contrast to the calculations in 1+1 di- 
mension (5; see also actual DLCQ calculations in 3+1 by 
Hiller |T2]. 

The reason is inherent to Dirac's relativistic vertex in- 
teraction (ki,hi\V\k2,h2',k3,fi3), in which some particle 
'1' is scattered into two particles '2' and '3' with their re- 
spective four-momenta k and helicities h, see Fig. ^ The 



matrix element for bremsstrahlung, for example, is pro- 
portional to k±, (ki,] \V\k 2 , |: fci, T) ' x |ki_|. scc Table 9 
in pQ, when the quark maintains its helicity while irradi- 
ating a gluon with four-momentum fcg = (xP + , kj_, k^). 
Singularities typically arise when over the squares of such 
matrix elements is integrated with repect to kL as in the 
integrations of perturbation theory. 

As usual, regularization is not unique, and can be done 
in many ways. But not all regularization techniques of the 
past are applicable in a Hamiltonian approach. Dimen- 
sional regularization, for example, is not applicable in a 
matrix approach which is stuck with the precisely 3+1 
dimensions of the physical world. Perturbative regulariza- 
tion |13) is not applicable in the non-perturbative context. 
The Fock space regularization of Lepage and Brodsky [H] , 
see also pQ, has blocked renormalization for many years. 
It was impossible to find suitable counter terms. In recent 
work, the invariant mass squared regularization has been 
abandoned in favor of Pauli- Villars regularization |12U14j . 
But thus far it is unclear how the approach |14| is appli- 
cable to the spectra of physical mesons. After applauding 
the light-cone approach [0)], Wilson and collaborators [05] 
have attempted to base their considerations almost en- 
tirely on a renormalization group analysis, but no specific 
technology has emerged thus far. 

In recent years I have come to favor another regula- 
rization, vertex regularization, which allows for renorma- 
lization and which is sufficiently simple on the technical 
level to be carried out explicitly. 

The singularities are avoided a priori by vertex re- 
gularization, by multiplying each (typically off-diagonal) 
matrix element with a regulating form factor F: 

(ki,h 1 \V\k 2 ,h 2 ;k 3l h 3 ) => 
(kuh^kzM-MM) F{Q) . (5) 

It took me several years to realize that it is the Feynman 
four-momentum transfer, Q 2 = —{ki — k 2 ) 2 , across a ver- 
tex, which governs an effective interaction. The minimal 
requirement for such a universal form factor is then 

F(Q) = F(Q;A) = {1> tQ 2 Z^. (8) 

It is satisfied by the step function, F(Q) = 0(Q 2 — A 2 ). 
The limit A — > suppresses the interaction all together, 
the limit A — > oo restores the interaction and its prob- 
lems. Any finite value of A 2 restricts Q 2 to be finite and 
eliminates the singularities. But the sharp cut-off gener- 
ates problems in another corner of the theory and F(Q) 
must be an analytic function of Q, as to be seen below. 

Vertex regularization takes thus care of the ultravio- 
let divergences. The (light-cone) infrared singularities are 
taken care of as usual by a kinematical gluon mass. 

Emphasized should be that vertex regularization does 
not chop off individual Fock space constituents, nor does 
it chop off Fock space components with large momenta, as 
in Fock-space regularization 0IH1 ■ Possible problems with 
gauge invariance are thus less severe and unlikely. Spoken 
in terms of matrices, the dimension of the Hamiltonian 



Hans-Christian Pauli: Solution to the non-perturbative renormalization of gauge theory 



3 



Fig. 1. 

Re- 

gu- 

la- 

ri- 

za- 

tion 

of 

the 

in- 

ter- 

ac- 

tion 

by 
ver- 
tex 
re- 
gu- 
la- 
ri- 
za- 
tion. 
In 
a 

ma- 
trix 
el- 
e- 
ment, 
as 
il- 
lus- 
trated 
on 
the 
left 
for 
a 

ver- 
tex, 
a 

quark 
changes 
its 
four- 
momentum 
from 

fci 
to 

fc2, 

i.e. 

Q 2 = 
-(fa- 
te) 2 - 

The 
ver- 
tex 
in- 
ter- 
ac- 
tion 
is 

reg- 
u- 
lated 
by 



4 



Hans- Christian Pauli: Solution to the non-perturbative renormalization of gauge theory 



matrix can be as large as it wants to be. The limit A 
is represented by a diagonal matrix. 



3 Renormalization of the total Hamiltonian 

The non-perturbative renormalization of the Hamiltonian 
was stuck for many years by the fact that the coupling con- 
stant g multiplies the regulator function F(Q) in Eq.0. 
It has always been clear that one may add non local coun- 
ter terms but is was not clear how they could be 
constructed. An explicit construction is given below. 

What follows is the general but abstract procedure of 
modern renormalization theory rT7j . 

Suppose to have solved the problem for a fixed value 
of the 7 'bare' parameters in the Lagrangian, the cou- 
pling constant g = g and the 6 flavor quark masses 
rrif — (to/)oj and for a fixed value of the exterior cut- 
off scale A = Aq. Suppose further that these 7+1 pa- 
rameters are determined such, that the calculated mass 
squared eigenvalues agree with the corresponding ex- 
perimental numbers. Next, suppose to change the cut-off 
by a small amount 6 A. Every calculated eigenvalue will 
then change by 6Mf. Renormalization theory is then the 
attempt to reformulate the Hamiltonian, such, that all 
changes SMf vanish identically. 

The fundamental renormalization group equation is 
therefore: 



dMf 



dMf 



g=g ,m f =mj ,A=A Q 







(7) 



for all eigenstates i. Equivalently one requires that the 
Hamiltonian is stationary with respect to small 5 A: 



5H, 



LC 



. 



(8) 



Hence forward reference to (go, rrif , Aq), to the 'renorma- 
lization point', will be suppressed. 

Next, I specify a non local counter term, as follows. To 
the original Hamiltonian one adds a counter term Hamil- 
tonian, whose interaction has exactly the same structure 
except that the regulator function F(Q) is replaced by an 
other function F(Q), that is by 



F(Q, A) F(Q, A) = F(Q, A) + C(Q, A) 



(9) 



Here, the regulating function F(Q) — F(Q, A) is supposed 
to be a given function. The counter term C(Q) — C(Q, A) 
is unknown, subject to the constraint that it vanishes at 
the renormalization point, 



C(Q,A) 



= 



(10) 



A=A 



The goal is then to find a C(Q,A) such that Eq.JSJ) is 
satisfied identically. 

According to Eq.© the light-cone Hamiltonian de- 
pends on A only through the function F(Q; A). Its varia- 
tion with respect to A is therefore 



SH 



LC 



5F- 



LC 



Eq.Q as the fundamental equation of renormalization 
theory can therefore be replaced by 



— dF 
SF = SA- = 0, 



(11) 



since the variational derivative of the Hamiltionian with 
respect to the regulator cannot vanish identically. 

The fundamental equation (|llfl defines then a differ- 
ential equation for C(A) — C(Q; A): 



dC(Q;A) _ dF(Q;A) 



(12) 



dA dA 
which, in its integral form, includes the initial condition 
A 



C(Q,A) = - Jds 



dF(Q,s) 
ds 



F(Q,Ao)-F(Q,A). (13) 



The renormalized regulator function is manifestly inde- 
pendent of A: 



F(Q,A) = F(Q,A)-{ 
= F(Q,A)1 
= F(Q,A ) 



C(Q,A) 
(F(Q,A Q )-F(Q,A)) 



(14) 



5F 



The counter term improved Hamiltonian is renormaliza- 
tion group invariant. By construction, the value of Ao must 
be determined by experiment. 

One should emphasize an important point: In deriving 
Eq. (|14|) , it was assumed implicitly that the regulator func- 
tion has well defined derivatives with respect to A. The 
thcta function of the sharp cut-off, however, is a distribu- 
tion with only ill defined derivatives; therefore, it must be 
removed from the class of admitted regulator functions. 

Another point is equally important: The step function 
has only a single parameter, A. Renormalization theory 
admits any function F(Q,A) = f(Q/A) as long as f(x) 
has well denned derivatives with respect to x — Q/A. 

For a long time, the main result of this section looked 
to me like sheer witchcraft: One selects a particular value 
of A = Ao, and then it turns out that the Hamiltonian 
is stationary with respect to small changes in A, for any 
value of Aq\ 

How can one believe this strange result? Is it only for- 
mal? — The abstract procedure remains somewhat aca- 
demic since a reliable numerical solution of the bound 
state equation, Eq.lJT}, would exhaust a large fraction of 
the worlds computer power, not to speak about the many 
variations which should be performed. To carry out such 
a programme would probably be impossible. Progress has 
come only from recent work on a particular model which 
did allow to formulate a paradigmatic example |18j in 
modern renormalization theory, see below. 

The result is in sharp contrast to perturbative renoma- 
lization, where the renormalized quantity is independent 
of the regularization scale. Here, in the non-perturbative 
renormalization of a Hamiltonian, the regularization scale 
Aq is fixed from the outset and becomes a member of the 
canonical parameters, like mass and coupling constant. It 
may not be driven to the extreme, to infinity or to zero, 
it must be kept at the finite value fixed ab initio. 
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4 The effective (light-cone) Hamiltonian 

In a field theory, one is confronted with a many-body prob- 
lem of the worst kind: Not even the particle number is con- 
served. In order to formulate effective Hamiltonians more 
systematically, a novel many-body technique had to be 
developed, the method of iterated resolvents 0E] , whose 
details are not important here. 

Important is that the method of iterated resolvents 
does not chop off Fock states or Fock state classes, and 
that it does not violate on purpose symmetries like gauge 
and Lorentz invariance. 

Important is as well that the effective light-cone Hamil- 
tonian -ff C Lc has the same eigenvalue as the full light-cone 
Hamiltonian H^c and that it generates the bound state 
wave function of valence quarks by an one-body integral 
equation in (x,kj_) [H]: 



i? m 2 . 



— y 

3tt 2 ^ 



x 1 — X 

dx'd 2 k\ 



faj y/x(l~x) X <(l-X<) 



x F(Q q )F(Q q 



a(Qq) , Oi(Qq) 

2Q 2 a 



2Q| 



x [^k 1 ,h 1 )^u(k' 1 ,h' 2 )][v(k 2 ,h 2 )^v(k 2 ,h 2 )] . (15) 

One has arrived at Eq.lTJJ: M 2 \^ qq ) = H^d^qq)- Here, 
M 2 is the eigenvalue of the invariant-mass squared. The 
associated eigenfunction iphihs ( x > Hi) is the probability 
amplitude (x,kj_,/ii;l — x, — kj_, h 2 \ 1 Fqq) for finding the 
quark with momentum fraction x, transversal momentum 
kj_ and helicity h\, and correspondingly the anti-quark. 
Expressions for the (effective) quark masses mi and rf%2 
and the (effective) coupling function a(Q) are given in [Hj. 
Q q and Q q are the Feynman momentum transfers of quark 
and anti-quark, respectively, and u(ki,hx) and v(k 2 ,h 2 ) 
are their Dirac spinors in the Lepage-Brodsky convention 
jS] , given explicitly in PP . They are arranged in the Lorenz 
scalar spinor matrix 

(/ii,/i 2 |sK,/i2> = Kfci^ihMfciX)] 

x[v(k 2 ,h' 2 )^v(k2,h 2 )} (16) 

which is a rather complicated (matrix) function of its six 
arguments x, x' , kj_, k^_, as tabulated in [TJJ. Finally, the 
form factors F(Q) restrict the range of integration and 
regulate the interaction. Note that the equation is fully 
relativistic and covariant. 

It should be emphasized that Eq. l|15fl is valid only for 
quark and anti-quark having different flavors ■ The 
additional annihilation term for identical flavors is omit- 
ted, but presently being investigated [20] ■ It should also 
be emphasized that the same structure was obtained by a 
completely different method, with Wegner's Hamiltonian 
flow equations In [2TJ it is also shown why the con- 
cept of a 'mean momentum transfer', Q 2 = \ (Q 2 + Q|) 
is a meaningful simplification. It allows to replace Eq. (|15|) 



by 



M 2 ip hl h 2 (x,k± 



+ 



rn 



X 1 



■>p hl h 2 (x,k.j_) 



rix'ri 2 k| VV i/4 (x',kl) q(Q) 

= 2~-K(yj 



3^ 2 fa J y/x(l-x)x'(l-x') 



x [5(ki I /n)A(*' lI 'i'i)]^,%''(*2.M] • (17) 

The form factors F(Q) have made their way into the reg- 
ulator function 



R(Q) = F 2 (Q) . 



(18) 



The same bilinear combination occurs also in the explicit 
expressions for to and a(Q) 6 , such that one can replace 
the F(Q)'s by R{Q) as an independent variable. 

It is a light cone pecularity that the (light-cone) Hamil- 
tonian -Hlc is additive in kinetic energy terms and the 
interaction pQ. The same holds for the effective Hamilto- 
nian H c g. Since it is allowed to subtract a c- number from 
i/ C Lc and to divide the result by another c-number, one 
can define an effective Hamiltonian H c g, implicitly by 



-HeLC = (mi + to 2 ) 2 + 2 (mi + to 2 ) H eS 
and explicitly by 



-Heff = 



H 



cLC 



(toi + m 2 y 



2 (toi + TO 2 ) 



(19) 



(20) 



Such a manipulation does not change the cigenfunctions. 
The eigenvalues then have the dimension of an energy 

#eff|*W) = E \&qq) , 

arriving this way at Eq.Q- As will be seen below, H c h 
has much in common with a conventional non-relativistic 
Hamiltonian in momentum representation, but the ap- 
proach remains fully relativistic and covariant. 

Note that mass and energy in the front form, on the 
light cone, are related by 



M 2 = (mi + m 2 ) + 2 (mi + m 2 ) E , (21) 



and not by 



M 2 = (mi + to 2 + E) 

= (mi + to 2 ) 2 + 2 (mi + m 2 ) E + E 2 
^ (mi+m 2 ) 2 + 2(m 1 +m 2 )E , (22) 

as usual in the instant form which quantizes the system 
at equal usual time. Only if the energy is negligible, i.e. 
only if (£/(mi +TO2)) 2 <C 1, the two relations coincide. 



4.1 The effective Hamiltonian in different clothes 

A rather drastic technical simplification is achieved by a 
transformation of the integration variable. One can substi- 
tute the integration variable x by the integration variable 
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k z , which, for all practical purposes, can be interpreted pQj 
as the z-component of a 3-momentum vector p = (k z , kj_). 
For equal masses mi = m 2 = m, the transformation is, 
together with its inverse, 



x(k z 



(23) 



(24) 



Inserting these substitutions into Eq. I|17|) and defining the 
reduced wave function iphih 2 by ^1 



^fo lft , 2 (a!,k]_) = — j== ' = <Phih 3 (h , kj_) , (25) 



(26) 



leads to an integral equation in 3-momentum space, in 
which all reference to light-cone variables has disappeared: 



For anti parallel helicitics hi = —h 2 (singlet) the model 
interpolates between two extremes: For small momentum 
transfer Q, the '2' in Eg. (|29l> is unimportant and the Cou- 
lomb aspects of the first term prevail. For large Q, the 
Coulomb aspects are unimportant and the hyperfine in- 
teraction is dominant. The '2' carries the singlet-triplet 
mass difference. 

For parallel helicities hi — hi (triplets) the model re- 
duces to the Coulomb kernel. The model over emphasizes 
many aspects but its simplicity has proven useful for quick 
and analytical calculations. Most importantly, the model 
allows to drop the helicity summations, simplifying the 
problem enormously on the technical level. 

The model can not be justified in the sense of an ap- 
proximation, but it emphasizes the point that the '2', or 
any other constant in the kernel of an integral equation, 
leads to singularities and thus to numerically undefined 
equations. Replacing the function a(Q) by the strong cou- 
pling constant a s — g 2 /Air completes the model assump- 
tions. 

For the singlet, the ST-model translates Ea. (|27|l iden- 
tically into 



Ef hl h 2 {p) = ~ <Phth 2 (p) 

2m r 



(27) 



, V K h- 2 (p') a c (Q)R(Q) 
y/A(p)A(p>) Q 2 4m 2 

x [u(ki,hi)ru(k'i,h' 2 )} [v(k' 2 , h' 2 )^v(k 2 ,h 2 )} 



5? £ h 

h',h' 



with a c (Q) = |a(Q). The reduced mass for mi = m.2 = m 
is m r = m/2. Its Fourier transform gives Eq.pJ. 

5 Two over-simplified models 

Krautgartner et al ^Oj and Trittmann et al rjTJ have 
shown how to solve an equation like Ea. (|17l) numerically 
with high precision. But since the numerical effort is con- 
siderable, it is reasonable to work first with an over sim- 
plified model, as specified next. 

The Singlet-Triplet model. Quarks are at relative 
rest when = and x — x, with x = TOi / (mi + m 2 ) . 
An inspection of Eq. (33) in JHj reveals that for very small 
deviations from the equilibrium values, the spinor matrix 
(hi, h2\S\h' 1 , h' 2 ) is proportional to the unit matrix, 

(hi,h 2 \S\h'ih' 2 ) ~ 4mim 2 S^j^ 6 h2th > 2 . (28) 
For very large deviations, particularly for kj 2 3> kj 2 , holds 
g 2 ~ki 2 , and (II \S\ Tl) c 2ki 2 . (29) 
The Singlet- Triplet (ST) model combines these aspects: 
(hi,h2\S\h' 1 ,h' 2 ) =Sh u h[ h 2 ,h' 2 (hi,h 2 \S\hi,h 2 ) ,(30) 

(31) 



(hiM\S\hiM) _ J ir ip+2, torhi = -h 2 , 



Q 2 



4miffl2 
Q' 2 ' 



for hi = h 2 



E V(P) = £~V(P) 



2m r 



(32) 



Oc_ 

2tt 2 



d 3 p' ( 4m 2 



R(Q) 

4m 2 



V(P') 



with The equation for the triplets is obtained by 

dropping the '2'. The first term in this equation, p 2 /2m r , 
coincides with the kinetic energy in a conventional non- 
relativistic Hamiltonian. This is remarkable in view of the 
fact that no approximation to this extent has been made. 
The fully relativistic and covariant light-cone approach 
has no relativistic corrections in the kinetic energy! 

Since the first term in Ea. l|32l) is a kinetic energy, the 
second must be a potential energy — in momentum rep- 
resentation. In principle, it could be Fourier transformed 
with c~ lpr to configuration space with the variable r. But 
due to the factor A(p)A(p') in the kernel, the resulting 
potential energy would be non-local. A Fourier transfor- 
mation with Q, i.e. with e~ lC * r , would not make much 
sense to me. 

To avoid this non-locality, I introduce the 'Fourier sim- 
plification' 



A( P ) EE 1 . 



(33) 



With A(p) = 1, the mean four momentum transfer Q 2 
reduces to the three momentum transfer q 2 = (p — p') 2 . 
In consequence, the kernel of Ea. H32|) . 



a {4m 2 t 2 \ R(q) 



2ir 2 \ q 2 



Am 2 



(34) 



depends only on q = p — p'. Its Fourier transform is a 
local function, 



V(r) = / d 3 q c- iqr U(q) 



(35) 
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playing the role of a conventional potential energy in the 
Fourier transform of Eq.J22l, i-e. in 



E ip(r) 



2m r 



V(t) 



(36) 



Here is the Schrodinger equation from Eq.Q in full glory! 
Despite its conventional structure it is a front form equa- 
tion, designed to calculate the light-cone wave function 

The jj-model. In an early stage of the present work, 
I thought that the Coulomb term in Eg. 13411 needs no re- 
gularization. In consequence I replaced the kernel by 



a ( 4m 2 



2R(q) 



Am 2 



(37) 



referred to as the T|- model. 

In view of the present insight, see [21], this misses 
the point completely and is completely false. Nevertheless, 
the model has played an important role in developing the 
present renormalization theory. 



6 Paradigms for renormalization theory 

In this section, I apply the tl-model to the the general 
procedure as outlined in Sect. [21 It is simple enough to 
go to every single aspect in concreto. In particular, I can 
solve the eigenvalue problem, H e Lc\<fi) = M 2 \ifi). 

The lowest eigenvalue Mq agrees with the pion mass 
at the given values of 



a = ctQ = 0.69 , 
m = mo = 1.16Z\ , 
A = A = 3.80Z\ , 



(38) 



with the mass unit A = 350 MeV. The triple will be called 
the renormalization point. 

When changing the parameters a and A (at a fixed 
value of m — mo), all eigenvalues M 2 will change and be 
functions of (a, A), among them the mass-squared eigen- 
value for the ground state, Mp (a, A). This function is dis- 
played as a contour plot in Fig. Oleft). Note that this 
plot is not schematic but the result of an actual calcu- 
lation. Calculated with 100 meshpoints in each direction, 
the figure represents some 10 solutions of the eigenvalue 
problem. It goes without saying that such can be done 
only for a sufficiently simple model, the tj-model, with a 
sufficiently simple soft regulator 



R(Q;A) = Ro(Q;A) 



A 2 



A 2 + Q 2 



(39) 



and with sufficiently fast numerical procedures. 

According to the general outline in Sec. [3| one must 
make sure that the invariant mass squared spectrum stays 
invariant, 8M 2 (a,A) — 0, for infinitesimal variations 8 A. 
I do that by construction, by introducing a new regulator 



R(Q,A) = R(Q,A)-(A 2 - A 2 ) 



dR(Q,A) 
dA 2 



(40) 



with the soft regulator R(Q, A) given by Ea. 1)39(1 . 

Both R and R have the same dimension. Both are in 
the class of admitted functions, see Eq.©. Both are equal 
at the renormalization point A = Aq: 



R(Q,A ) = R(Q,A ) 



(41) 



The derivative of R vanishes at the renormalization point 
by construction, since 



dR(Q,A) 
dA 2 



2Q< 



(A 2 



Al) 



(A 2 + Q 2 f 



(42) 



This derivative is certainly not the zero function for all A, 
and renormalization can be valid at most locally, in the 
vicinity of the renormalization point. 

Based on the Hellmann-Feynman theorem, one expects 
a vanishing derivative of the eigenvalue with respect to A 
near the renormalization point. The actual numerical re- 
sults in Fig. fright) illustrate this very convincingly: All 
contours have zero derivatives at A — Aq. The Hamilto- 
nian is locally renormalized: 6Mf(a,A) ~ 0. 



6.1 Renormalization by a sequence of regulators 

With Rq(Q, A) given by Ea. l|3l?|l . the above can be put into 
a systematic form. I can define a whole series of regulator 
functions R n (Q,A) by 



Ri(Q,A) 
R 2 (Q,A) 



-Al) 



Ro(Q,A) - Q 2 ^ i+ Qiyi , 



Ri(Q,A)-Q 



(A 2 +Q 2 ) 3 



(43) 



Rn(Q,A) = Rn-xiQ,*) 



q 2 _r^gr 



(A 2 



Each member is a genuine regulator satifying the minimal 
requirements of Eq,©. All are different from the original 
regulator in Ea. H39l) . and none can be preferred over the 
other. Their derivatives given by 



dR x {Q,A) 

_dA 2 
dR 2 (Q,A) 
dA 1 



dR n (Q,A) 
dA 2 



OQ2 (A 2 -A 2 ) 

— ^ (y 1 2 + Q2)3 

- on 2 ( a . 2 - a q) 2 



(n + 1)Q 



2 (A 2 -A 2 r 



they all are stationary at the renormalization point. Car- 
rying on the procedure to higher and higher orders, the 
contours corresponding to Fig. ffi right) become broader 
and broader. In the limit of large order the contours be- 
come flat: One has realized the fundamental renormaliza- 
tion group equation: dM 2 (a,A) = 0. 

One can show this even analytically by considering the 
limit n — > oo. With 



A 2 



A 2 



A 2 + Q 2 ' 



(44) 
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the regulators R n (Q,A) behave in the limit n — > oo like 
R n (Q, A) = R (Q, A) - Q 2 4 ~ 5o [l + x + x 2 + 



A 2 



A 2 + Q 2 
A 2 - A 2 1 



A 2 + Q 2 A 2 + Q 2 l-x 
A 2 A 2 A 2 



A 2 + Qi A 2 + Q2 A^ + Qi 



A 2 



A 2 + Q 2 



(45) 



which is manifestly independent of A. This result repre- 
sents a beautiful and pedagogic example for how the re- 
normalization group works. 



6.2 Renormalization by counter terms 

The regularization by counter terms has been discussed 
in Sec. I repeat it here for clarity. Introduce a counter 
term Hamiltonian by 

R(Q, A) => R(Q, A) = R(Q, A) + C{Q, A) . (46) 

The counter term C(Q, A) is the unknown, subject to the 
constraint 



C(Q,A) 



. 



A = An 



The renomalization condition, 

d~R{Q,A) 
dA 



0, 



(47) 



(48) 



(49) 



defines the derivative of C: 

dC{Q:A) _ dR(Q;A) 
dA ~~ dA ■ 

Including the constraint by the integral form, 

A 

C(Q,A) = - fds dR ® ,s) = R(Q,A )-R(Q,A) , (50) 



take the Fourier transform according to Eq. l|35|l . and de- 
fine this way precisely the ill-defined Hamiltonian of the 
'Coulomb plus delta' interaction according to Eq.JSSJ), i.e. 



H 



2m r 



a 
r 



5^(1 



(53) 



Such a problem could then be solved by his 'subtraction 
method,' see f.e. Frederico et al [22] , a method based on 
a T-matrix approach like in scattering theory. 

Frederico's initiative has put us into "the unique po- 
sition to compare two drastically different schemes, both 
conceptually and numerically, and show that they agree" 
(citation from [22]). 

This is of course no proof, but: 'A reasonable man is 
satisfied with evidence, only a fool needs proof,' as the 
famous mathematician Marc Kac has said in the context 
of trying to prove Boltzman's H-theorem. 

The step to the 'oscillator plus delta' interaction [22|, 



H = ^-a+ f -r 2 +^(v) 



(54) 



was then straightforward and gave beautiful agreement 
with all of the the experimental meson spectra. 



7 Perturbative versus non-perturbative 
renormalization 

Perturbative renormalization addresses to render quanti- 
ties of physical interest independent of the regularization 
scale. Conceptually, the scale can then be driven to infin- 
ity in the remainder of the theory, A oo, such that it 
drops out from the theory. Perturbative renormalization 
is adequate for scatterings problems. 

Non-perturbative renormalization says: 'No, don't do 
that. Do not drive A to infinity.' Set the universal renor- 
malization scale Aq in accord with experiment and give 
it the same rank as the other physical parameters of the 
theory, the coupling constant a and the mass m. Non- 
perturbative renormalization is adequate for bound state 
problems. 



determines the counter term. The renormalized regulator 
R(Q,A) = R(Q,A ) (51) 
is then manifestly independent of A, in analogy to Eo. 1)14(1 . 

6.3 Renormalization by the subtraction method 

After my talk at one of the Trento meetings, Tobias Fred- 
erico approached me with the proposal to put R(q) = 1 in 
Kq.OJ- 



. , a f 4m 2 \ 1 



27r 2 V q 2 



Am 2 ' 



(52) 



8 Summary and Conclusions 

This work is an important mile stone on the long way from 
the canonical Lagrangian for quantum chromo dynamics 
down to the composition of physical hadrons in terms of 
their constituting quarks and gluons, by the eigenfunc- 
tions of a Hamiltonian. 

As part of a on-going effort, a denumerable number of 
simplifying assumptions had to be phrased for getting a 
manageable formalism 3 . Among them is the formulation 
of an effective interaction by the method of iterated resol- 
vents [B]. The strongest assumption in the present work, 
the simplifying Singlet-Triplet model in Sec. ^ is also the 
easiest to relax, as shown in the next paper. 
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The biggest progress of the present work can be found 
in Sects. [3 and It is related to a consistent regulariza- 
tion and renormalization of a gauge theory. The ultraviolet 
divergences in gauge theory are caused less by the possi- 
bly large momenta of the constituent particles, but by the 
large momentum transfers in the interaction. In a Hamil- 
tonian approach, such as the present, one has not much 
choice other than to chop them off by a regulating form 
factor in the elementary vertex interaction. 

The form factor makes its way into a regulator func- 
tion which suppresses the large momentum transfers in the 
Fourier transform of the Coulomb interaction. The arbi- 
trariness in chopping off the large momentum transfers 
is reflected in the arbitrariness of the potential at small 
short distances. The arbitrariness allows for a pocket in 
the potential which binds the quarks in a hadron |24| . 

The problem is then how to fix this function with its 
many parameters, by experiment. In practice this is less 
difficult than anticipated. Actually, it has ben done al- 
ready It suffices to determine only three parameters, 
two continuous ones and one counting index. 

The present work opens a broad avenue of applications, 
among them also the baryons and physical nuclei. But 
much work must be done in the future before such a simple 
approach be taken serious. 
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